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ex. 

Abstract: We compute the parity violating part of the time-dependent gravitational 
response function of an ideal gas of Weyl fermions up to third order in the derivative 
expansion and give its full tensorial structure. Our main results are two functions that 
parametrize the energy-momentum tensor in terms of gauge-invariant combinations of vec- 
tor and tensor metric perturbations. The zero frequency limit of these functions is related 
with the anomalous constitutive relations and with the full anomalous partition function in 
the presence of gauge and mixed anomalies. In particular, our results imply the existence of 
a previously unknown invariant contribution to the parity-odd partition function at third 
derivative order that we explicitly construct. Beyond the static limit, the gravitational 
response function may provide valuable insights into time-dependent phenomena driven by 
anomalies. 
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1 Introduction 

One of the most interesting and surprising developments in thermal field theory has been 
the recent realization that anomalies induce unexpected modifications in the constitutive 
relations of relativistic hydrodynamics [1-3]. As a consequence, there are allowed transport 
phenomena [4-6] related only to anomalies which manifest themselves by parity violation. 
An important feature of the new anomalous terms is that they do not break time-reversal 
symmetry and hence the equilibrium response to a perturbation is not accompanied by 
an increase in entropy. In this regard, it has been argued that a systematic derivation of 
constraints on the constitutive relations should be possible solely from the knowledge of a 
consistent functional describing the equilibrium thermodynamics, without having to make 
use of an entropy current [7-9]. More specifically, this functional has been constructed 
for an anomalous charged fluid, to second order in the derivative expansion in 1+1 dimen- 
sions [10], and to third order in 3+1 dimensions [11]. These constructions, supplemented by 
the appropriate Bardeen-Zumino terms [11-13] to be added to the consistent currents, pro- 
vide higher-order anomalous contributions to the constitutive relations due to gravitational 
and mixed anomalies. 
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As mentioned above, an important feature of the anomalous contributions to the con- 
stitutive relations is that they break invariance under parity. To see the implications of 
this fact, let us consider the non-relativistic limit of the anomalous constitutive relation 
connecting the vorticity to the momentum density [14, 15] 

9 = XpvVxv, (1.1) 

where v is the fluid velocity in a generic frame (not in the Landau frame). According to the 
theory of hydrodynamic fluctuations [16], the non-zero static susceptibility xpv is given 
by the following sum rule 

XPyochm ^ — / <ko— , (1-2) 

in terms of the imaginary part Xggi^t Q) °f the retarded momentum- momentum correlator. 
For a translationally invariant system this correlator takes the general form 

X%(",<1) = q ^TXL^,q) + - ^) X T(u,q)+ie ijk q k XPv(uJ,q). (1.3) 

Since all the momentum density components g l have the same signature under discrete 
symmetries, the last term of (1.3) can be different from zero only if the equilibrium state 
used in the definition of x l gg{t,x,0,O) = ( [g l (t, x), g J (0, 0)] ) is not invariant under parity. 
For a system of chiral fermions this requires a nonzero axial chemical potential, and xpv 
is expected to be an odd function of this potential. This observation suggests that the 
study of the parity violating piece of the two-point function for gravitons may be useful 
in order to deepen our understanding of anomalies in hydrodynamics and also to gain 
information about time-dependent phenomena driven by anomalies, a field that remains 
relatively unexplored. 

In this paper we undertake the computation of the parity violating part of the gravita- 
tional response function of an ideal gas of Weyl fermions at small frequencies and momenta 
compared to the chemical potential and temperature. In other words, we compute the grav- 
itational response function in the hard-dense-loop approximation. The procedure is based 
on previous work by Rebhan [17], who has given the full tensorial structure of the leading 
corrections oc T 4 of the graviton self-energy. A brief review of the work by Rebhan, gen- 
eralized to non-vanishing chemical potential, is presented in the next section, where the 
resulting energy-momentum tensor is shown to correspond to the constitutive relations of a 
perfect conformal fluid at lowest order in the derivative expansion of the fluid fields. This is 
extended to third order in the derivative expansion in section 3, where our main results are 
two functions that parametrize the energy-momentum tensor in terms of gauge-invariant 
combinations of vector and tensor metric perturbations. The zero frequency limit of these 
functions is related with the anomalous constitutive relations in section 4, and with the 
form of the full anomalous partition function in the presence of gauge and mixed anomalies 
in section 5, where we find a previously unknown invariant contribution to the parity-odd 
partition function. Our conclusions and possible applications of this work to the study of 
time-dependent hydrodynamic phenomena driven by anomalies are presented in section 6. 
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2 Leading contribution at low momentum 



In this section we briefly review previous work by Rebhan [17] on the graviton self-energy 
at leading order, and connect it with the constitutive relations of a perfect conformal fluid. 
The graviton polarization tensor is defined by 



W vpa {x-y) = 



where T is the effective action and 



9=T) 



~g{T p °{y))) 



(2.1) 



(T pu ) = ^_/ r . (2.2) 



-g <ty 



/"' 



Note that eq. (2.1) implies 

5 (V^T^Oc))) =-\\ #yW<»{x - y)h pa (y), (2.3) 

where the retarded version of Jl p,upa '(x — y) has to be used in order to compute the corre- 
sponding induced change to linear order in h pu = g pu — 7] pu . 

It was shown in [17] that, for small momenta \q \,q <C Tq,/iq the leading behavior of 
n^" pa (Q) can be written in terms of the integral 

r d 3 k K P K V K p K a 

I^iQ) =T J2 J (2?r )3 K 2f K + Q) 2 . = + Q p = (iv n ,q), (2.4) 

where v n (uj n ) are bosonic (fermionic) Matsubara frequencies, and the pure vacuum (Tq = 
//o = 0) divergence has been subtracted. Concretely, the following combination of indices, 

W pa (Q) = 2I^ pa {Q) - r] a p (l a ^Pr] ua + I^P^ + I a ^°r, up + l a ^rf p ) , (2.5) 



produces the correct graviton polarization tensor satisfying the Ward identities implied 
by general covariance. After the sum over Matsubara frequencies is performed, the small 
momentum behavior of I puXcr (Q) can be extracted by making the rescaling q° = iv n — > 
eq°,q — > eq, and keeping only O(e ) terms. The integrand obtained in this way is pro- 
portional to the energy density which, for a left (or right)-handed Weyl field, is given 

by 

e = ^k(n F (k - mo) + n F (k + Mo)) = ^ + ^ + A (2.6) 

The remaining angular integration yields the entire dependence of IP uXa on q°/\q\- In 
particular, the explicit computation of I pupcr (Q) in the static limit (q° = 0) yields the 
following non-vanishing components [17] for the polarization tensor 

n 0000 = -5e, n 00ii = - e , n ° 404 = --, 

3 

11"'-' - n ii00 = -e, nijij = £ (2.7) 

3 3 
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where the value of IP'JJ is valid only for i ^ j, with YL nu = e/Yl. 

An important feature of these results is that its form exactly corresponds to the con- 
stitutive relation of a perfect conformal fluid at lowest order in the derivative expansion of 
the fluid fields 

(T^) = (e + P)u»u v + Pg^. (2.8) 

This can be seen by considering the system in a curved background. The most general 
static metric which is preserved by the Killing vector dt may be written as [9] 

ds 2 = - e 2<7{s) (dt + aj {x)dx j ) 2 + g ij (x)dx i dx j . (2.9) 

Now, in a comoving coordinate system, the fluid velocity becomes = <5ge -cr , and the 
temperature and chemical potential also acquire a dependence on the position given by [18] 

T(x) = H?oor 1/2 T = e- CT T , 
H{x) = e no- 

Here To and no are constants which may be viewed as the temperature and chemical 
potential in the absence of gravity 1 . Thus, in the conformal case P = e/3, the induced 
corrections to (T^) are simply obtained by making the replacements e — > e~ 4r7 e, u v — > 
e _<J 5o and jf' v —> g^ u . To linear order in this yields 

5(- 
H 

H 
H 

where /loo = — 2<r, Ziofc = — o-k and hij = gij — 5ij. A look at (2.7) shows that the coefficients 
in (2.11) exactly agree with the components of H^ upa in the static limit. Thus, they encode 
the form of the constitutive relations for the energy-momentum tensor at the lowest order 
in a derivative expansion of the fluid fields. In the remaining of this paper we will extend 
Rebhan's computation to higher orders in the momenta and analyze the implications for 
the constitutive relations and parity violating partition function. 



g(T™)) = ^h 00 + b -Y / h kk , 



-g(T 01 )) = 3%, 
=£<T*» = £ -h m 



(2.11) 



IP*)) 



3 Parity-odd response function 

As explained in the Introduction, our interest will be in the parity violating part of the 
gravitational response function. From the fact that (T pv ) in flat space-time is of the 
form (2.8), which is parity-even, it follows that, in order to obtain the parity-odd part of 
5(T pu (x)), eq. (2.3) can be simplified to 

5(T^(x)) = - l -j An^(* - y)h pa (y), (3.1) 

1 Note that To and /io are the same constants that appear in the partition function Z = Tr e~ t - H ~ t± ° N ^ T ° 
even if the Hamiltonian includes the static gravitational field. 
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where the required response function is given by 

). (3.2) 



n^^( x -y) = -i 6(x° - y°) < [T^ \x\T"° (y)] ) - 2 ( i(V 9^)^ (x)) 



g=T)l 



For an ideal gas of left-handed Weyl fermions, the first term in (3.2) takes the following 
form in the imaginary time formalism 

/J 3 k 

1 7Fj (3.3) 
x W°{K + Q, K)] K2{K \ Q)2 , K° = lUJn + (J,, 

where V- = (1 — 7s)/2 and the fermion-fermion-graviton three- vertex, which can be read 
from the energy- momentum tensor T^ v in Minkowski space (A. 3), is given by 2 

V» V {K,P) = - 4 [ 1 »{K + Py + ^{K + Py] - + (3.4) 

Up to parity even contributions (see appendix A for details) the second term in (3.2), 
coming from the seagull diagram, can be written 



i " r d 3 k i (3 - 5) 



where a va = j[7 I/ ,7 <r ]. Here the zero subscript from T and fi has been omitted. 

3.1 Leading contribution to the parity-odd response function 

As reviewed in the previous section, for momenta <C \fJ-\,T, the leading behavior of 

rP t/pcr is proportional to the energy density and can be written in terms the prototype in- 
tegral (2.4). The next to leading order contribution to Ii^ upa in the momentum expansion, 
which is odd in the chemical potential and linear in Q, turns out to be governed by the 
integral 

I^(Q)-TT f ^ K " KURX (3 6) 

1 (Q) - 1 ^ J (27r)3 R2{K + Q)2 • (3.6) 

Concretely, the leading parity-odd contribution to 111 is obtained by picking the 75 
in the projector V— = (1 — 7s)/2 together with as many Ks as possible in the numerator 
of (3.3). This gives 

nf P V, q) = \QcJ){i\ (e a ^I vrTX + e^Pl^ +{p^a). (3.7) 



2 The gamma matrices obey {7", 7"} = -2^, 75 = -^^poYYYl" where e i23 = -e 0123 = -1. 
Note that we use the mostly positive signature. 
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Similarly, the parity-odd part of H2 is given by 



where 



nf pff ((? u ,q) = -^Qavpx (e^o^H* + e^V*^) + (P (3-8) 



r d 3 k K p . . 

"" = r E/p^- <«■<» 

As before, the evaluation of (3.6) involves Matsubara summation and rescaling of the 
resulting integrand through q° — > eq°, q — > eq. The integral obtained by keeping the O(e ) 
term, which we will denote lj^ X , has all components proportional to 

f°° dk 1 
/ (0)=/ ^k 2 (n F (k-fi)-n F (k + f i)) = —(^ + 7T 2 fiT 2 ). (3.10) 

" 

In particular, we have 

T} afl lffi = H^ = -I {Q) u\ (3.11) 



where u 7 = (1,0) is the velocity of the fluid in the local rest frame. Using this relation 
to eliminate H x in favor of and adding the contributions of IIi and H2, gives the 
following expression for the leading parity-odd contribution to the polarization tensor 

njyv,*/) = sQavpx****^ - l^^f) + ... 

+ e «f>>>P(S%-r ] >»rh«)+(P 

where 6% = 8»8% + 8*6%. 

We can check this result by noting that the covariant conservation law for the energy- 
momentum tensor imposes constraints on the tensorial structure of U pup(T (Q). Since the 
covariant derivative involves a combination of derivatives of 5{T pu ) and Christoffel symbols, 
and both of them are of the same order in h pa , it follows that the vanishing of 8{T pu (x)) 
gives rise to Ward identities connecting Q fl H fJ-upcr (Q) with the one-point function (T a @) [19]. 
The consequences of this requirement in relation with (2.4) were fully analyzed in [17]. With 
regard to the parity- violating part of H^ upa ^Q), the conservation law imposes uniquely 
Q ( ,n f "' ' pa (Q) = 0, since the parity-violating part of (T a ^) vanishes. Similarly, conformal 
invariance requires t] lJll ,Il p ' 1 ' pa (Q) = for this part. Now, U^ pa as given by (3.12) is 
obviously traceless, while transversality follows from the following property satisfied by the 
componets of I puX at leading order 

ij^Qx = \(Q^ + Q»8»)r) a pl$\ (3.13) 

which is easily proved from the results in table 2 of appendix B. Actually, up to an overall 
normalization, the structure of the leading parity-odd response function (3.12) is uniquely 
determined by the conditions of transversality and tracelessness together with (3.13). 
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In order to evaluate II^ pcr we need explicit expressions for all the components of I^ uX . 
This is conveniently done by using the decomposition I^ uX = ^ ■ cfT^ vX in terms of the 



six symmetric tensors 



fivX 

1 

fiuX 



fiuX 



rpjlvX 
rpfluX 

1 5 

rpfluX 



Q^Q U Q\ 

QVQ u u x + Q»u v Q x + u»Q u Q\ 
Q»u v u x + u' 1 Q u u x + u' 1 u u Q x , 

u»u u u x , 



(3.14) 



The contractions of I^ uX with the basis tensors are given in table 2 of appendix B. Inverting 
these relations gives the coefficients Cj 



ci(q°,q)/I(o 
C2(q°,q)/I(p 
c 3 (q°,q)/I {0 
C4(q°,q)/I(o 

c5(q°,q)/i(o 

C(j(q°,q)/I(o 



5q° 3g°(5Q 2 ~ V) r , 

V 4? ' 

5Q 2 -2g 2 3Q 2 (5Q 2 -4 9 2 ) 



4g 4 4g 6 
g °(5Q 2 + 2g 2 ) 15gV Q 



£(<A</), 



Aq 2 + V (? ' 



(3.15) 



Q 2 



2 /5Q 2 + 2g 2 15Q 



Aq 4 



+ 



Aq e 



-L(q°,q)), 



Q 2 7-0 + 



3Q' 



4g 2 4g 4 



£(<7°,<?) 



where L(q°,q) = Qi(q°/q) is the Legendre function of the second kind, which results from 
using the retarded prescription iv n + + — > q° in the analytic continuation of the integral 



t 



dt- 

! iv n - qt 



-L(q°,q) = - 

q q 



-1 + 2- In 
2q 



q° + q 



'i 9 -m 

2 q 



0\2 



(3.16) 



This result can be greatly simplified by using the constraints imposed by (3.13), which 
restrict the number of independent functions to three, for example c\,C2 and C3. Clearly 
ci does not contribute, and we are left with only two combinations of the functions C2 
and C3. Due to the transversality of the response function, the result can be written in a 
particularly transparent way in terms of the two projectors Pf and Py, 



jf J,U 



V 



(u-Q) 2 + Q 2 



[u ■ Q(u"Q u + u u Q v ) + Q"Q" - Q 2 u»u v ] 



Q 2 



(3.17) 



T • 
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Then the response function (3.12) adopts the simple form 



+ icj(q°, q)u a Q p [e a ^P^ + e a ^P^ + (p O a)] , 
where the functions cl,t are given by 

q 2 

cv(q°,q) = ^o(9° c 2 + c 3 ) 



^+*V*) (l + ^L(g°,g)), 



24vr 2 
1 
2 

96vr 2 



CT( ( ?°,g) = ^[Q 2 C2-<? c 3 ] 



(3.19) 



In the static limit these become 

<*(o,ff) = -^ 8 + *Vi*), (3 . 20) 
or(0,g) =0. 

These are the main results in this subsection. 

3.2 Parity-odd response function at higher orders in the momenta 

In order to use the response function as a source of constraints on the anomalous constitu- 
tive relations at higher order in the derivative expansion, we will need H^ pa (Q) up to third 
order in Q. The form of the second order contribution comes from the explicit contribution 
of order Q 2 to the trace in (3.3) together with the piece of 0(Q) in I a Pl . From the explicit 
formula (3.5) for II2, it is obvious that the seagull diagram does not contribute beyond first 
order. Therefore the second order correction reads 

where we have introduced the integral 

d 3 k K a KP 



I^(Q)-TT I Sk ™" (3 22) 



which in this case is needed only at O(Q ). A computation using the results for the 
integrals in appendices B and C shows that, although neither I?£ nor I 1 ^ 1 vanish, the 
specific combination in (3.21) does. Therefore the second order contribution to the parity 
odd gravitational response function vanishes identically. 
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The third order contribution may be expressed similarly. Besides at 0{Q 2 ) and 
I a @(Q) at 0(Q), we also need the O(Q ) piece of the integral i^, which turns out to be 
independent of the temperature 



k" 



d 3 k 



(2tt) 3 K 2 (K + Q) 2 



8tt 2 



91 

r,2 



(3.23) 



The appropriate combination is now 



(3.24) 



Using the explicit results for the integrals in appendices B and C, and following the method 
outlined in the previous subsection, we find that the third order response function is still 
given by eq. (3.18), with the coefficients replaced by 



1927T 2 
,2 



Of 



/" 



1927T 2 

In the static limit they become 



2Q 2 3QV-2Q 2 ] 
—5- + 1 HQ » ?) 



Q 4 3Q 6 o 
V+2^'^. 



(3.25) 



cy(0,g) 
ct(0,<?) 



P 2 
2^ ' 



192vr^ 

V- 2 



(3.26) 



This completes the computation of the parity-odd response function to third order in Q v . 



4 Energy-momentum tensor and metric perturbations 

In this section we will use the results obtained for the parity-odd response function to 
derive the general form of the parity violating part of the energy-momentum tensor. We 
will devote special attention to the static case, where our results can be compared with 
recent proposals in the literature. As shown above, the effects of metric perturbations on 
the energy-momentum tensor 



6{T 



in^(<A<z)V> 



(4.1) 



can be parametrized by the two independent functions cy(q°,q) and ct(<?°,<?) 



ip^V, q) = iCv ( g °, q) {u .^2 + Q 2 u «Qp V a ^ p PV + f^'P? + {p^ *)] (4 2) 



■l(q°,q)u a Qp [e^'if" + e a ^P^ + (p O a)} . 
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Scalar 


Vector 


Tensor 




-2a 






hoi 


-dib 


(SO 






c&ij + didjd 


diFj + djFi 


h^ 



Table 1. SO(3) components of a general perturbation of the metric, where we have defined a\ 
and a\ = dib as the solenoidal and irrotational parts of Sgoi = —ai(t,x). 



It should be noted that, as a consequence of the form of IP 1 ' po ', the response of the 
energy-momentum tensor depends only on gauge-invariant combinations of the metric dis- 
turbances. In order to simplify the analysis, a general metric perturbation has been decom- 
posed into (SO (3) irreducible components in table 1. Note that the vector fields a\ and 
Fi are solenoidal, while h^ is traceless and satisfies dih^ = 0. Direct substitution shows 
that scalar perturbations do not produce any parity- violating effect on T pu . For vector 
perturbations, not necessarily time-independent, the change in the expectation value of 
(T^ v ) depends only on the combination m + dtFi, which is gauge-invariant 

5(T *) = cy(q ,q)ie^ k q 3 (-a k +iq°F k ), 

Sip*) = cy (<A q) iq° (e mn q m q J + e jmn <T<f) (-a» + iq°F n ), 

where g J = q 3 /q. Thus, cy(q°,q) parametrizes the response to vector perturbations of the 
metric. Similarly, cj(q°,q) parametrizes the response to tensor perturbations hij, which 
are gauge-invariant by construction 

6(1**) = -ar(q ,q)e ilm 6l n iq l h mn + (t o j). (4.4) 

These are the main results in this paper. In what follows, we will compare them with other 
results in the literature. 

4.1 Static limit of response functions and anomalous constitutive relations 

Arguments based on linearized hydrodynamics [16] show that the small velocity field of the 
fluid v 3 (t, x) is the quantity that plays the role of external force coupled to the momentum 
density l° l in the perturbing Hamiltonian 

*- = -/*»*"<*.•*(*.»>■ <«) 

Hence a comparison with 

= , (4.6) 

enables us to identify hoi = — with the fluid velocity v l in the static limit q° = 0. Thus, 
using (3.20) we find the following constitutive relations in the static limit at leading order 
in the momenta 

8(T (H )=c v (0,qyJ k iqiv k , 
5(T ij ) =0. 
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We see that cy(0, g) may be interpreted as a parity-violating susceptibility connecting the 
vorticity with the momentum density. In fact, in hydrodynamics with quantum anomalies, 
the momentum density to linear order in v [14] is precisely given by 

9 =(£|P)n»Vxi), (4.8) 

when one uses a frame where the entropy current does not have anomalous part i.e., where 
Jg = su^. This agrees with our result (4.7) which, by (3.20), implies the following value 
for the anomalous susceptibility 

XV = c V (0, q) = 0" 3 + ^V^) ■ ( 4 - 9 ) 

Actually, the value of the parity violating susceptibility has been related to the gauge and 
mixed anomaly coefficients in the anomalous conservation equations 

V„ J£,v = \^ vp " {Zc A F^F pa + c m R a ^RP ap a) , (4.10) 
1 

— ( 

2 



V V TZ = F\r cm + \c m V u {e^F pa R^ afi ) , (4.11) 



through the relations [11] 

Xv = 2(c M (iT 2 - CAfJ 3 ), c 4d = -%ir 2 c m . (4.12) 
The values of the anomaly coefficients for a left-handed spinor in (3 + 1) dimensions, 

c A = 8c m = — ?-2, (4.13) 

then imply our value (4.9) for the anomalous susceptibility. 

Similarly, we may use (3.26) to find the following 0(Q 3 ) corrections to the constitutive 
relations in the static limit 



S(T oi ) = J^2^%a k , (4.14) 
5t ^ j ) = -i£-2<? msik V 2 dih km + (t O j). (4.15) 

1927T Z 



Note that (4.14) simply gives a correction to the momentum density (4.8) 

S 9 = "TT^ v2 ( V x «) = -c m /"V 2 (V x v). (4.16) 

1927T Z 

On the other hand, (4.15) describes a qualitatively new effect. For a tensorial perturbation 
depending only on z 

fh + (z) h x (z) ()\ 
h x (z) -h + (z) , 

y o o 0y 



(4-17) 



eq. (4.15) gives 

T 11 = —T 22 = 2cmntiZ(z 



(4.18) 
T = -2c m fih + (z). 
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Predictions for the effects of gauge and mixed anomalies on the constitutive relations 
at higher orders order in the momenta have only recently been given in (1 + 1) [10] and 
(3 + 1) dimensions [11]. In order to connect our results with those in [11] we collect some 
of their formulae. Using a frame where the energy-momentum tensor takes the form 

(TH = ( £ + P)v?u v + PgT + u ^Qa + u v Qa + tT, (4-19) 

with u^q^ = u fl T^ u = 7^ =0, the following results were obtained [11] at 0(Q 3 ) in (3+1) 
dimensions 

<t A = -2fic m v% = -2iic m e^ pa u v V ' p R aa u a + ... (4.20) 
= ±iicj% = Aiic m & a< »e v> ^u p V a R a p + . . . (4.21) 

where the dots stand for contributions that vanish at linear order in the metric perturba- 
tions and are thus invisible in our computation. In this formula 

V <^> = AWA W 7 W - ^-A pv A pa V fKJ (4.22) 

3 

and A^ = g pv + u^u v '. In order to compare with (4.14), we evaluate the pseudovector v£ 
at linear order in the metric perturbation. This yields v® = and 

4 = ~\^ k {djV 2 h k0 - dodjdihtk) + 0(h 2 ), (4.23) 

where we have taken u 11 = 5q + 0(h). Upon substitution in this expression of the metric 
components in table 1 we obtain 

q\ = -2fic m 4 = —^2^ 2d Mk + d t F k ) + 0(h 2 ), (4.24) 

which differs in sign from our contribution (4.14) to 5(T° l ) in the static limit. On the 
other hand, note that precisely the combination at + dtFk, which we have argued should 
be identified with the fluid velocity, appears in this expression. Proceeding similarly with 
the pseudotensor gives the non-vanishing components 

f = ~-/ m {did a d a hj - d0d a h m a ) + (t O j) + 0{h 2 ). (4.25) 

Taking the static limit and substituting the metric components in table 1 gives 

r l l = ^cJi = --J^e^dthJ + (i o j), (4.26) 

which exactly agrees with (4.15). In the next section we will explore the relations of our 
linearized results with the full anomalous partition function given in [11] and will be able 
to explain the sign discrepancy noted above. 
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5 The anomalous partition function and the static response function 



In this section we consider in more detail how our previous results for the static response 
functions fit with known facts about gauge and gravitational anomalies of the underlying 
theory, particularly with the form of the full anomalous partition function. As we will see, 
an important feature of the third order response (3.26) is that it reveals the existence of an 
invariant contribution to the partition function which, in principle, does not seem to arise 
from general arguments about anomalies. 

We will assume that an external time-independent gauge field A = Ao(x)dt + Ai(x)dx l 
is present besides the static gravitational field given in (2.9). In this background with 
Killing vector = (1,0), the equilibrium temperature and chemical potential become 
position-dependent. They are defined in terms of the invariant length of the time circle 
and the Polyakov loop Pa as [7] 

•l/To 



Jo J o 



M(*0 r „ „ f 1/T ° , ,_„,u M*) 



In 



P A = / dTA^W 

Jo 



T{x) J ^ ' T ' 

where Tq 1 is the length of the compactified imaginary time. 
5.1 The first order generating functional 

Let us start with the first order in the derivative expansion. Using the definition of the 
energy-momentum tensor in terms of the partition function 3 

the linearized constitutive relation (4.7) may be rewritten in terms of Ao and To (or equiv- 
alently \x and T) 

5(T 0t ) = -T 5 -^- = {2c A Al ~ 2c M T 2 AoW jk d J a k + ..., (5.3) 

where e 123 = 1, and e ljfc = A/53 will denote the corresponding tensor. Since the 
quantities dj and A% transform like the components of a covariant vector under spatial 
diffeomorphisms and Aq behaves like a scalar, it is clear that W must include the terms 

- y J d 3 x^e ijk A 3 aidjak + c 4d T J d 3 x^e ijk A aidja k . (5.4) 

To reconstruct the full dependence of W on the gauge field, we complete the action with the 
requirements of gauge invariance up to a U(l) anomaly, and invariance under Kaluza-Klein 
transformations [9]. These correspond to redefinitions of time, t — > t' = t + <f>(x), without 
change in the spatial coordinate, and preserve the form of the metric if ai transforms 



3 Here W = InZ — —Q/To, where Z is the grand partition function and f2 refers to the thermodynamic 
potential, a functional of the background quantities depending solely on x. 
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as 5ai = —di4>. Under such a transformation the gauge field changes as 5Aq = 0, SAi = 
—CfpAi = —Aoditj), and the combination Ai — aiAo remains invariant. As this kind of gauge 
invariance may be viewed as a manifestation of the underlying diffeomorphism invariance, 
which is not anomalous at the first derivative order, it is natural to impose this requirement 
on the partition function. A short computation shows that, up to a total derivative, the 
resulting Kaluza-Klein invariant action is given by 



W = -^ [ d 3 x^e^ k A 3 a i d j a k + ^ [ d 3 'yfax e^ k , J&A i d j a k 
J-o J -toy 

~ljf / d3x V93 eijkA oA i d j A k -c M T J d 3 x^e ijk (Ai - A ai)d jak . 



(5.5) 



Thus the static gravitational response, together with the requirements of invariance under 
Kaluza-Klein transformations and gauge invariance up to a £7(1) anomaly, determine the 
partition function at first order in the derivative expansion. The first three terms in (5.5) 
constitute the anomalous part of the consistent partition function, while the last one is 
gauge invariant. Our result (5.5) agrees with the form of the parity-odd partition function 
proposed in [9] for a CPT invariant theory. The consistent gauge anomaly follows from the 
variation of this action under a gauge transformation 5Ai = d{A, 5Ao = 0. This induces a 
change 

5W = ^-j cFxy/toAeWdiAodjAk, (5.6) 

which shows that the consistent anomaly is precisely determined by the cubic term in the 
static momentum density correlator. The variation of the consistent current under such a 
gauge transformation is given by 

5J° = -2c A -^LdjA k diA, 



8J { = -^A-^LdjAodkA. 

Therefore the anomalous gauge invariant current is obtained by the compensating shift 

J c ^ ov = ° + 2c A e^ pa A u d p A a . (5.8) 
V~9 oA^ 

Note that the Bardeen-Zumino term in this equation contributes to the current at first 
order in the derivative expansion. 

5.2 The third order generating functional 

Now we turn our attention to the connection between the static gravitational response at 
third order in the derivative expansion and the mixed anomaly. First, we introduce the 
Pontryagin density defined by 



V = \*R^R v ^p, (5.9) 
where the dual Riemann tensor is given by 

*R» u a ? = L a ^ T R^ pT . (5.10) 
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Locally, the Pontryagin density can be written as a total divergence 



VpK* = V, (5.11) 
where is the Chern-Simons topological current 

Ra = e^T^d^ + |l^rk) • (5.12) 



We have seen from (4.15) that the response to a tensor perturbation of the form (4.17) 
is given by 



T ll = _ T 22 = 2Cm ^ {z) 

T 12 = -2c m ^h%{z). 



(5.13) 



This response involving only spatial indices may actually be derived from the quadratic 
portion of the functional 4 

W K = ~ f d 3 xV=9~AuK>*, (5.14) 



To 

which, for Ao = A 4 , A% = 0, becomes 

^quad _ Cm 



W' K - 

J-o 



J d 3 x n(h' + (z)h^(z) - hl(z)h' x (z)). (5.15) 

But Wjf Uad , when evaluated for the vector perturbation hoi = —&i{z) , ho2 = —ai(z), 
ho3 = 0, 

wF* = - c fJ = ~WoS ^rW'Ki*) - °i(*K(*)). (5-i6) 

does not produce the required response. Indeed its variation yields 

SWT** 



T 01 = -To^- = c mf ia^(z), 

(5.17) 



6a% 



i-TT^quad 

rr02 rji ovv K in < \ 

T u = -T = -CmHa-L (z), 

which is minus our result (4.14). This is consistent with our observations at the end of 
last section. Therefore, an additional contribution to the action will be needed in order 
to properly account for the vector response. The following is a natural, possibly unique 
choice that preserves U{1) gauge invariance, three-dimensional diffeomorphism invariance 
and Kaluza-Klein gauge invariance 

W inv = ^ [ d 3 x^e^ k g m ^d m {^f mn )d 1 {A k - a k A ), (5.18) 



4 In the linear approximation the required components read T 11 = ToSW^** / 'dh+ = — T 22 and T 12 = 
T <5^, uad /^x. 
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where the inverse metric g kl is used to raise the lower indices of = diOj — djOi, and 
ci is a constant to be determined shortly. In terms of differential forms the integrand is 
proportional to d(A — A$a) A 5da, where 5 denotes the codifferential. The virtue of W- mv 
is that the linearized tensor response remains unaffected, while its contribution the vector 
response reads 

5W {2) 

-To— ™- = -2c lM 4"(*), 

5a \ 2) (5-19) 
-To-r^- = 2dK '(z). 

00,2 

If we could identify Wk with the anomalous partition function, then c\ = c m would 
be the appropriate choice in order to reproduce our linearized results (4.14). However, Wk 
can not be identified with the partition function. Instead, the non-invariant functional Wk 
plays the role of a local counterterm that interpolates between two alternative definitions 
of the consistent effective action. These two choices preserve either diffeomorphism or 
gauge invariance. Following Bilal's notation [13], if we denote by I^ 1 ) the diffeormorphism 
invariant effective action for the mixed gauge- gravitational anomaly, and by the gauge 
invariant one, the interpolation is given by 

W K + = T (2) . (5.20) 

We must then consider the contributions to the energy-momentum tensor from I^ 1 ' 2 ) and 
the relation with our results for the static response. 

The complete anomalous contribution to the partition function at third derivative 
order, together with its implications on the constitutive relations, have been recently sta- 
blished in [11]. With the vorticity and acceleration of the fluid given by 



the authors of [11] construct the well-behaved covariant current 

1 



(5.21) 



j£ = -4W^ pa u u u p L0 a - + IRpvUPxf - 2a x a x - ^ x lo x J uf, (5.22) 

where W^ vpa is the Weyl tensor and R pa is the Ricci tensor. This current, evaluated for 
the equilibrium fluid velocity = e~ a V^ in the background (2.9), has the same divergence 
as the topological current 

= V. (5.23) 

This implies that the local functionals 

J f (5.24) 

v{2) = %j d * x ^Mf m -Kn, 
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satisfy (5.20) and have the required properties: is obviously invariant under three- 
dimensional diffeomorphisms, but not gauge invariant, whereas, thanks to (5.23), is 
gauge invariant under 5Ai = <9jA, 5Aq = 0, but not diffeomorphism invariant. As a 
consequence, these functionals give rise to two different types of consistent observables. 
Let's first consider the combination + W- mv . It is easy to see that the component 
is cubic in the metric perturbation and does not contribute to the linearized T^ v . Thus, 
for Aq constant and A% = 0, the quadratic part of matches that of Wk, and has 
identical contributions to the linearized T^ u . On the other hand, the Bardeen-Zumino 
term required to compensate for the lack of invariance under diffeomorphisms vanishes if 
the gauge field strength T a p does [11], and the linear response results will agree with the 
variation of the quadratic part of + W{ nv , which coincides with that of Wk + W{ nv , for 
the choice c\ = c m . 

For a different choice of the generating functional, namely for I^ 1 ) + W\ nv , the con- 
clusions are identical, although they require a little more work. This functional may be 
viewed as the third derivative counterpart of the first derivative action (5.5), since both of 
them are diffeomorhism invariant. Differentiation with respect to A^ yields a gauge and 
diffeomorfism covariant anomalous current with no need for a compensating term. But, 
while the first derivative term SW/Sg^ from (5.5) was gauge invariant, the consistent con- 
tribution ^(rW + W- my )/Sg l j il , is not, and a Bardeen-Zumino tensor T^ z has to be added 
in order to produce a gauge invariant and generally covariant energy-momentum tensor 



2T <5(r«+W- inv ) 



T»» = , 7 + Tg z . (5.25) 



As 

Jm i s cubic in the metric perturbation, rW does not contribute to the linearized 
T^ u . Thus the third-order derivative response, which obeys the Ward identity of general 
covariance, must precisely match the linear portion of the Bardeen-Zumino tensor, together 
with the contribution from Wi nv . This turns out to be the case if c\ takes the value c\ = c m , 
which is determined from (5.19). Indeed, an explicit computation shows that the results 
in (4.14) and (4.15) may be rewritten as 

c)W {2) 

6(Tn = 2T — ^ + 2c m d 3 (XrVo + XrVo) • (5.26) 

The non-linear generalization of the second term, 

Tg z = 2c m V A (*R x ^ up A p + *R Xu ^A p ) , (5.27) 

exactly matches the form of the Bardeen-Zumino energy tensor given in [11] when the 
anomalous piece of the generating functional is precisely I^ 1 ) . 

Thus the results at the third derivative order obtained in this paper show that a 
consistent generating functional describing the effects of the mixed anomaly must include 
an additional invariant piece W- mv given by (5.18) with c\ = c m , a feature that ultimately 
can be traced to the opposite signs for cy and cj in (3.26). 
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6 Discussion and outlook 



In this paper we have computed the subleading corrections to the time-dependent parity 
violating graviton response function up to third order in the derivative expansion. In 
particular, we have obtained the complete invariant decomposition of the response function 
on the appropriate tensor basis (3.14) and (C.l). We have then exploited the transversality 
of the response function to rewrite our results in terms of the two invariant functions cv,t 
which condense all the relevant information. 

From the static limit of the response function we have extracted the anomalous con- 
stitutive relations at first (4.7), and third order (4.14), (4.15) in the derivative expansion, 
which we have compared with recent results in the literature. We have also shown that 
one can reconstruct the complete parity violating partition function at first order in the 
derivative expansion (5.5) from the linearized corrections to the energy-momentum ten- 
sor (4.7). The situation at third order in the derivative expansion is more involved, but 
we have shown that our expressions for the constitutive relations are fully compatible with 
very recent results on the form of the anomalous partition function at that order if one 
includes the previously unknown parity-odd invariant contribution given by (5.18). 

We have also obtained the complete dependence of the quantities cy,T on q°/q, from 
which one can compute chiral effects on the time-dependent departures from equilibrium 
8{T^ u (t, q)). For the parity-even part of the time-dependent response to metric perturba- 
tions, an analysis including a comparison with the results from the Boltzmann equation 
has been given in [20]. Now, having at our disposal a set of results for the parity-odd 
part of the response, we can pose the question about the relationship between the field 
theory approach presented in this paper and a possible kinetic description involving the 
background metric. In this regard, the modification of the non-equilibrium kinetic equa- 
tion that take into account chiral magnetic and anomalous Hall effects has been recently 
obtained in [21, 22] in the presence of an external electromagnetic field. In addition to this, 
it would be interesting to take advantage of the results at the linear level in this paper to 
obtain a Vlasov-type equation from which one could derive non-equilibrium chiral effects 
caused by metric perturbations or weak curvature backgrounds. We leave the consideration 
of this issues for future work 5 . 
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A Expansion of the action and seagull terms 

In this appendix we obtain eq. (3.5) for the seagull contribution to the parity-odd part 
of the response function. The action for a fermion in a curved background is given by 

5 After this work was completed, we became aware of [23], which deals with issues similar to the ones 
mentioned in this paragraph. 
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S = f cfixyJ—gC where 



c = i [^v^ - (v^jyv] , (A.i) 



with V^ip = d^tp — Tfiij}. The spin connection is related to the vierbein e y a by 

r M = lh a n b Ke bU]fl = ^7 6 ]<0V^ - r^)e/?«, (A.2) 



where greek and latin letters are used for curved and Minkowski indices respectively, with 
{7^,7"} = —2g^ v and {7 a ,7 b } = —2r] ab . Expanding the action S = J d 4 xC in powers of 
hfu> = g^v — 1J/J.U gives S = So + S\ + S2 + ■ ■ ■ where So is the action in flat space-time, 
Si = \ J d^xT^h^ with 



T -- 
L w ~ 4 



_ ^ y _ ^ y ~\ % — ^ y 

V'Tm d u ip + VHV d ^ - - r q^uipl a d a ip, (A. 3) 



and 



S 2 = — J d 4 x^ {a^,Y} *l> V aP K^ p hp u 

~\f d A x (3^T^ - 2^T^) h ap h pv (A.4) 

+ i J d\ {ifPrT - ^yj^a^^. 

with a^ CT = jpy 1 ', 7°]- The three- vertex fermion-fermion-graviton in eq. (3.4) then follows 
from the Fourier transform of eq. (A. 3). Similarly, the seagull contribution in eq. (3.5) 
follows from eq. (A.4) through 

+ <A5) 

o 

- \W\T UV ) +rfP(T a ^)), 

where we have used the fact that the the equations of motion imply (jC) = 0. As (T^ u ) in 
flat space-time takes the form (2.8), which obviously preserves parity, this establishes the 
validity of eq. (3.5). 
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QaQpQ-ylffi /I(n 
U a 7lfrlffi/I {n 




n 



1 



3g°Q 4 
-3g°Q 2 
2q 2 -3Q 2 -2Q 2 L(q°,< 
4Q 2 

3g° + 6g L( (? ,g) 
-4g° 



n = 2 


-6Q 4 L(g°,g) 
12 g °Q 2 L(g°,g) 


2g 2 - 24(g°) 2 L(g°, 




Table 2. Contractions of -f^ 7 with the tensors in eq. (3.14) for n = 0, 1, 2. 



B The integrals P"' 

Table 2 gives the contractions of I pup with the basis tensors in eq. (3.14). The constants 
I( n ) are related to the Fermi distribution as follows 





poo 


dk 2/ 


1(0) 


Jo 


4^ k ( 






00 dk 






4T 2 






f°° dk 


1(2) 


= 48 J 


o 4^ 



1 



12^ 2 



(/i 3 +7rVT 2 ), 



1 Z" 00 djfc 1 
4 / T^ /c ( n ^( /c -/ 1 ) + n ^( fc + / x )) = 96^2^ 3 ^ 2 + 71-2712 



n F (/c - n) - n F {k + /x)) = ' 



(B.l) 



where the prefactors have been chosen for convenience. These relations can be inverted to 
give the coefficients in the expansion I^ vp = ^ ■ CjT pup ', with the following results: 



The tensor expansion of I pup at leading order has been given in eqs. (3.15) and (3.16). 
Tensor expansion of I pup at 0(Q) 



Cl(q°,q)/I { 
C2(q°,q)/I( 
C3(q°,q)/I( 

C4(<7°, <?)//( 

c 5 (q°,q)/I( 
c 6 (q°,q)/I( 



3 3(3Q 2 -2q 2 ) ^ 
^2 + W.?). 

2g° 6g°Q 2 o 
2 — L <? ,?), 

Q 2 3q°Q 4 (5Q 2 - 2q 2 ; 
q 2 q & 

_^ 2 

q 



L(q°,q), 



2 L(q°,q), 



0. 
0. 



(B.2) 
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Tensor expansion of I^ vp at 0(Q 2 ) 
A Cl (q°,q)/I {2 



c 2 (g°,g)// ( 

G4(?°,?)/I( 
c 5 (q°,q)/I {2 
C6(q°,q)/I(2 
C The integrals /' 



q°(5Q 2 -2 g 2 ) 3g° (5Q 4 - W + 8 g 4 ) r . . 
s H s M? ,9), 



9" 9" 
Q 2 (5Q 2 -4g 2 ) 3Q 2 (5Q 4 -6g 2 Q 2 + V) 
g 4 + g6 

5g°Q 4 3g°Q 4 (5Q 2 - 2g 2 



L(q°,< 



q°Q 2 W r , o 
~^ ^ L((? 

5Q 6 3Q 6 (5Q 2 - 2g 2 ) 

Q^_ 3Q 1 
^2 



L(g c 



Table 3 gives the contractions of with the following basis tensors 

rf = w, 

it = *r, 

with the constants Ii n ) given in eq. (B.l). 





n = 


n = 1 


QaQpI^J I{n+l) 


-2Q 2 







2q° 


12Q 2 


u a u p I*/I {n+1) 


2 + AL(q°,q) 


-24q L(q°,q) 


T a P 1 T 

Vapl {n) /l{n+l) 


-4 






(B.3) 



(C.l) 



Table 3. Contractions of with the tensors in eq. (C.l) for n = 0, 1. 



The coefficients of the expansion P IU = ^ ■ cfT^ v are given by 

• Tensor expansion of at leading order 

o 2 6Q 2 -4g 2 r , o 
ci(q ,q)/I(i) = —2 ~i HQ >?)> 

C2(<A ?)//(!) = ^ + ^rW> 9), 
C4(q°,q)/I {1) =^L( q ,q). 
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• Tensor expansion of at 0(Q) 



ci{q ,q)/I(2) 

C2(q°,q)/I(2) 
C3(q°,q)/i(2) 

Ci{q°,q)/I{2) 



0. 



0. 



24q L(q°,q), 
12Q 2 L(q°, q), 



(C.3) 
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